Sisyphus cooling and amplification by a superconducting qubit 
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Laser cooling of the atomic motion paved the way for remarkable achievements in the fields of 
quantum optics and atomic physics, including Bose-Einstein condensation and the trapping of atoms 
in optical lattices. More recently superconducting qubits were shown to act as artificial two-level 
atoms, displaying Rabi oscillations, Ramsey fringes, and further quantum effects^'^'^. Coupling such 
qubits to resonators*'^'^'^ brought the superconducting circuits into the realm of quantum electro- 
dynamics (circuit QED). It opened the perspective to use superconducting qubits as micro-coolers 
or to create a population inversion in the qubit to induce lasing behavior of the resonatoi^'^iiSiii. 
Furthering these analogies between quantum optical and superconducting systems we demonstrate 
here Sisyphus coolingi^ of a low frequency LC oscillator coupled to a near-resonantly driven super- 
conducting qubit. In the quantum optics setup the mechanical degrees of freedom of an atom are 
cooled by laser driving the atom's electronic degrees of freedom. Here the roles of the two degrees 
of freedom are played by the LC circuit and the qubit's levels, respectively. We also demonstrate 
the counterpart of the Sisyphus cooling, namely Sisyphus amplification. 

For red-detuned high-frequency driving of the qubit the low-frequency LC circuit performs work 
in the forward and backward part of the oscillation cycle, always pushing the qubit up in energy, 
similar to Sisyphus who always had to roll a stone uphill. The oscillation cycle is completed with a 
relaxation process, when the work performed by the oscillator together with a quantum of energy 
of the high-frequency driving is released by the qubit to the environment via spontaneous emission. 
For blue-detuning the same mechanism creates excitations in the LC circuit with a tendency towards 
lasing and the characteristic line-width narrowing. In this regime "lucky Sisyphus" always rolls the 
stone downhill. Parallel to the experimental demonstration we analyze the system theoretically and 
find quantitative agreement, which supports the interpretation and allows us to estimate system 
parameters. 



The system considered is shown in the inset of Fig. 1. It consists of a three-junction flux qubit-'^^, with the two qubit 




FIG. 1: (a) The energy levels of the qubit as a function of the energy bias of the qubit £{fx) ~ 2^oIpfx. The sinusoidal current 
in the tank coil, indicated by the wavy line, drives the bias of the qubit. The starting point of the cooling (heating) cycles 
is denoted by blue (red) dots. The resonant excitation of the qubit due to the high-frequency driving, characterized by Qro, 
is indicated by two green arrows and by the Lorentzian depicting the width of this resonance. The relaxation of the qubit is 
denoted by the black dashed arrows. The inset shows a schematic of the qubit coupled to an LC circuit. The high frequency 
driving is provided by an on-chip microwave antenna, (b) SEM picture of the superconducting fiux qubit prepared by shadow 
evaporation technique. 
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states corresponding to persistent currents of amplitude Ip flowing clockwise and counterclockwise. When operated 
in the vicinity of the degeneracy point, = <i>x/$o — 1/2 ~ 0, where $x is the magnetic flux applied to the qubit 
loop and $o = h/2e the flux quantum, the Hamiltonian of the qubit in the basis of the persistent current states reads 

H = -]^eUx)<yz-\^<y. . (1) 

Here Gx^cTz are Pauli matrices, A is the tunneling amplitude, and s{fx) = 'i^olpfx is the energy bias. The energy 
levels of the isolated qubit, separated by AE{fx) — \/£'^{fx) + are shown in Fig. 1. The qubit is driven by a 
high-frequency field with frequency and coupled via a mutual inductance M to a low-frequency tank circuit with 
frequency wt — 27r x 20 MHz <C A/h. Both circuit to the qubit can be included at this stage via their contributions 
to the external flux $2,. Since the eigenfrequency of the tank circuit, ujt, is much lower than the level spacing of the 
qubit, the oscillations of the current in the tank circuit can be treated in an adiabatic approximation, i.e., the current 
in the tank circuit shifts the bias flux of the qubit by ^x{t) = M[Idc + Iri{t)]. 

The system mimics the Sisyphus mechanism of damping (cooling) and ampliflcation (heating) known from quantum 
opticsi^. This mechanism is illustrated in Fig. 1. Here we describe the damping (cooling, hence marked in blue) for 
a situation where the driving is red-detuned, hcod < AE; the amplification (marked in red) for blue-detuning can be 
described in an analogous way. The oscillations of the current in the tank circuit, Iif{t), lead to oscillations of e{fx) 
around a value determined by the dc component, Idc- In the first part of the cycle, when the qubit is in the ground 
state, the current shifts the qubit towards the resonance, AE = hujd, i.e., the energy of the qubit grows due to work 
done by the LC circuit. Once the system reaches the vicinity of the resonance point, the qubit can get excited, the 
energy being provided by the high-frequency driving field. With parameters adjusted such that this happens at the 
turning point of the oscillating trajectory, the qubit in the excited state is now shifted by the current away from the 
resonance, such that the qubits energy continues to grow. Again the work has to be provided by the LC circuit. The 
cycle is completed by a relaxation process which takes the qubit back to the ground state. The maximum effect is 
achieved when the driving frequency and relaxation rate are of the same order of magnitude. If the relaxation is 
too slow the state is merely shifted back and forth adiabatically during many periods of oscillations. Note that the 
complete cycle resembles the ideal Otto-engine thermodynamic cycle^'*. 

Two types of measurements have been performed on the system. In the first the LC tank circuit is additionally 
driven near-resonantly by a low- frequency rf current, and the response of the LC circuit is detected using lock-in 
techniques. In these measurements we identify the infiuence of the high-frequency driven qubit on the effective quality 
factor and eigenfrequency of the tank circuit. We associate a reduction of the effective quality factor with cooling and 
identify regions in parameter space where the effect is optimized. 

In the second type of measurements the low-frequency r/-driving is switched off, while the emission of the LC 
tank circuit is monitored by a spectrum analyzer. This analysis probes the influence of the qubit on the effective 
quality factor, and, most importantly, allows us to determine the energy stored in the tank circuit, i.e., its effective 
temperature. Thus we are able to demonstrate cooling or heating of the LC oscillator. 

We start with a qualitative analysis of the flrst type of measurement. The LC tank circuit is driven by a current 
source with amplitude I^. Then the amplitude, Vr, of the voltage oscillations across the tank circuit, which is 
measured in the experiment, is 

Vt = LJj'Lx-^T Q (^T^T-^T ' (^) 

where It = ^2{I'^^) = QI^ is the actual amplitude of the rf current in the inductance, and Q is the effective quality 
factor of the tank circuit. It is given by the ratio, Q = 2ttWt/A, between the energy stored in the tank, Wt = LtIx/'^^ 
and energy loss per period A. The latter is composed of two contributions, A — At + Agis- The intrinsic losses of the 
tank are given by At = 2'kWt /Qo, where Qo is the intrinsic quality factor of the tank circuit. The average work done 
by the tank on the qubit in one period, Asis, can be estimated as follows: We consider the optimal situation when the 
oscillator brings the qubit into the resonance at the turning point of its trajectory (Fig. 1). We assume further that 
the state of the qubit is instantaneously "thermalized" , with equal probabilities to remain in the ground state or to 
get excited. In the latter case, after leaving the area of resonance the qubit can relax with rate Ft?,. The probability of 
qubit relaxing during the time < < < T (here T — 271 /lot) within an interval dt is given by dP — exp(— Ffli) V ndt. 
Thus we obtain the average value of the work (taking into account the probability 1 /2 of the initial state to be excited) 

Asis = MIjJt j dP [I - cos{ujTt)] 

= MIplTf{TR,UT) , (3) 



where 

/(r„..ri.(i-.— «'-)-^^ (4) 
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FIG. 2: (a) Amplitude of the rf tank voltage as a function of the dc magnetic bias of the qubit, f^'^, for a microwave driving 
frequency ujd = 2n x 14.125 GHz and various amplitudes of the rf current driving the tank circuit, (b) The height of the peaks 
(solid squares) and dips (solid circles) as a function of the voltage amplitude of the unloaded tank circuit Vto- The voltage Vro 
is equal to, e.g., Vr taken at /^'^ = 0.02 where the effect of the interaction between the tank circuit and the qubit is negligible. 
The height saturates near 200 nV. 



This function demonstrates that the optimal situation for damping or amphfication is reached when Tji ^ ujt- For 
the effective quality factor we obtain 

where ± stand for the damping/amplification. Substituting into Eq. ^ we arrive at 

Vt - Vto = TMcjTlpQofirR, ujt)/^ , (6) 

where Vro = QqujtLtIt is the voltage on the tank circuit far from the resonance. Eq. ^ provides the estimate for 
the increase/decrease of the tank voltage for large driving currents I^, such that the qubit spends most of the time 
away from the resonant excitation area. For weaker driving, when the times spend away and within the excitation 
area are comparable, the damping effect becomes weaker. 

Our results of the first type of experiment are shown in Fig. 2. The dips (peaks) correspond to the Sisyphus 
damping (amplification) of the tank circuit, i.e., to the decrease (increase) of the effective quality factor Q. The 
central dip is due to the shift of the oscillator frequency when the qubit is at it's degeneracy point^-^ and is not related 
to the Sisyphus effect. We observe from Fig. 2b that the additional voltage saturates for large driving amplitudes at 
approximately 200 nV. Using this value and Eq. ^ we see that indeed the relaxation time is close to the period of 
oscillations, 1/Tr w 0.95r. 

The results for the second type of experiment are shown in Fig. 3. The comparison shows that at the bias point 
corresponding to maximum Sisyphus damping the resonant line widens in good agreement with the results for the 
driven LC circuit. At the bias point corresponding to maximum amplification the line narrows, showing a tendency 
towards lasing behavior. We can extract the quality factors in these two regimes, as well as for the bias point far from 
the resonance where no additional damping occurs. The comparison of the experimental results with Eq. ([5]) yields a 
quantitatively reasonable agreement. 

Finally, the results displayed in Fig. 3 allow us to estimate the efficiency of the cooling or heating in the two regimes. 
Integrating the power spectra we observe that in the damping regime the number of photons in the LC circuit decreases 
by about 8% compared to the undamped case. The cooling effect was observed also for other microwave frequencies 
and lies in the range of 6% — 13%. Thus there is only little cooling. One reason is that our system is optimized 
towards maximum Sisyphus damping rather than minimum cooling temperature. Indeed, damping is optimized when 
the resonant point is reached at the turning point of the oscillators trajectory. In this regime even a small reduction 
of the oscillators amplitude, i.e., a weak cooling, suffices to switch the whole mechanism off. We expect to be able to 
optimize in future experiments the system towards efficient cooling as well. 

To develop the discussion further we model the Sisyphus damping as an additional effective bath coupled to the 
oscillator, which we characterize by a quality factor Qsis and temperature Tgis. The standard analysis gives the total 
quality factor and temperature, = Qq^ + Qgj^ and TQ~^ = TqQq^ + TsisQgis, respectively. Here Tq is the 
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FIG. 3: (a) Spectral density of the voltage noise in the LC circuit, Sv{i^), measured when the low-frequency r/-driving is 
switched off while the high-frequency qubit driving is fixed at uid = x 6 GHz. Sv is shown as a function of frequency (vertical 
axis) and normalized magnetic flux in the qubit f!^'^ (horizontal ajcis). (b) The integral, J Sv{'-o)dui, evaluated using the data 
from panel (a) as a function of f^'^ . The integral is directly proportional to the number of quanta (effective temperature) of the 
tank circuit. At optimal dc bias f^'^ the effective temperature T of the oscillator is lower than the temperature To away from 
the resonance, (To — T)/To = 0.08, corresponding to an 8% cooling, (c) Spectral density Svi^^) measured at three different 
values of /^'^ corresponding to damping (blue), amplification (red), and away from the resonance (black). These values of f!^'^ 
are marked, respectively, by blue, red, and black arrows in panel (a). 



temperature of the oscillator without Sisyphus damping, which in our case is determined mostly by the amplifier. 
Our data show that Q ^ Qo/2. Thus the coupling strength of the Sisyphus mechanism is comparable to that of the 
rest of the environment, i.e., Qgis ~ Qa ~ 340. In the present experiment it looks as if Tgis were comparable with Tq 
and we have only little cooling. 

In contrast, in the amplification mode of our experiment the line gets much sharper, which formally can be expressed 
by choosing Qsis negative, Qsis ~ —600. This means that if no other bath was present the system would show the 
lasing instability. Consistently, the Sisyphus mechanism, acting as an active medium with population inversion, is 
described by a negative temperature Tsis. Its value is not known, but we can get a lower bound for the effective 
temperature from the inequality T/Q > Tq/Qq. It predicts an increase in temperature and number of photons by 
a factor ~ 2. The experiment shows a 36% increase, which demonstrates the expected trend, even though some 
quantitative discrepancy remains. 

We now outline the theoretical analysis of the problem. Quantizing the oscillations in the tank circuit we arrive at 
the Hamiltonian 

+ huJTCL'a + gaz{a + a^) ■ (7) 

where the coupling constant is g = MIplT.o, and It,q = y^TiuiT /2Lt is the amplitude of the vacuum fluctuation 
of the current in the LC oscillator. The third term of Eq. ([7|) describes the high-frequency driving with amplitude 
Hrq. After transformations to the eigenbasis of the qubit and some approximations appropriate for the considered 
situations (and described in the Appendix) the Hamiltonian reduces to 

H = AE az + hQfio cos (uJdt) cos(^ + hcjT i'' a 

9 

Q 2 

+ gsm(az {a + a^) - —cos'^ Cf^z {a + a^) , (8) 



with AE = + and tanC = e/A. Thus we obtain the effective Rabi frequency ^IjiqCos( and the effective 
constant of qubit-oscillator linear coupling gsinC. As we need both these terms for the Sisyphus cooling the qubit 
should be biased neither at the symmetry point nor very far from it. The second order term oc g^ is responsible, e.g., 
for the qubit-dependent shift of the oscillator frequency^^ . 

To account for the effects of dissipation we consider the Liouville equation for the density operator of the system 
including the two relevant damping terms. 



P = -T [H,p] +Lqp + LrP 



(9) 
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FIG. 4: Color lines: experimental data. Black lines: numerical solution of Eqs. © for Qro — 271 x 0.5 GHz, Fr — 1.0 • 10*s ^ , 
r* = S.O-lO^s'V The central dip is due to the quadratic coupling term in ((8} which causes a shift of the oscillator frequencjJ^. 



As far as the qubit is concerned we consider spontaneous emission with rate Tji and pure dephasing with rate F* . 
Hence, we have 

+ -f {^zP^z - P) ■ (10) 

We neglect the excitation rate since the qubit 's energy spUtting exceeds the temperature. Thus the standard longi- 
tudinal relaxation time Ti is given by T^^ = F/j. The rates Tn and F* may depend on the working point, i.e., on 
^, which in turn is determined by the driving frequency ojd by the condition hujd ~ AE. The rates should describe 
the dissipation around these values of (. It should be mentioned that pure dephasing is frequently caused by the 1/f 
noise^, for which the Markovian description is not applicable. We hope, however, that the main features are still 
captured provided F* is chosen properly. 

The resonator damping term can be written in usual formii, 

LrP= {Nth + 1) {2apa'^ - a) ap - pa'' a) 

+—Nth(2a^pa~aa^p~paa'^, (11) 

where k — lot/Qo characterizes the strength of the resonator damping, and A^th — [G:X'p{hujT / ksT) — 1] ^ is the 
thermal average number of photons in the resonator. We solve the master equation ^ numerically in the quasi- 
classical limit, i.e., for n = (a^a) ^ 1, and for various choices of the unknown parameters F^^, F* , and riuo- As the 
system is harmonically driven, we determine the response of the observables/density matrix at the driving frequency, 
and, finally, find the amplitude of the driven voltage oscillations across the tank circuit. As shown in Fig. 4 by fitting 
the system parameters within a reasonable range we reproduce well the experimental findings. We should mention 
that a good fit can be obtained for a relatively wide range of values of the Rabi frequency Qrq and pure dephasing 
rate F* , provided we keep their product roughly constant. On the other hand, F/j is determined rather accurately. 

The Sisyphus damping by a superconducting qubit analyzed in this paper can also be applied to nanomechanical 
systemsiSiii^'^Si^. With proper optimization towards minimum cooling temperature one could cool a mechanical 
resonator which would enable, e.g., ultra-sensitive detection applications like MRFM^^ with quantum limited precision. 
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Appendix 

After transformation to the eigenbasis of the qubit the Hamiltonian ([7]) becomes 

H = AE az ~ TiflRQ cos (ujdt) {sinC (Jz - cos (fJx) 
+ hujTo'a + g {sm( az — cos( ax) {a + a^) , 

with tan ( = e/A and AE = ^/e'^ + A'^. Because of the large difference of the energy scales between the qubit and the 
oscillator, AE ^ Tlcut, we can drop within the usual rotating wave approximation (RWA) the longitudinal driving term 
—?i57flo cos (wtfi) sin C cr^. On the other hand, we retain the transverse coupling term, — gcosCtT^; (a + a^), but trans- 
form it by employing a Schrieffer- Wolff transformation. Us — exp (iS), with generator S — [g/ AE) cosC (a + a^) Cy, 
into a second order longitudinal coupling. The Hamiltonian then reduces to the form 



